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A model of photo-detection using a Bose-Einstein condensate in an atom-chip based micro trap is analyzed. 
Atoms absorb photons from the incident light field, receive part of the photon momentum and leave the trap 
potential. Upon counting of escaped atoms within predetermined time intervals, the photon statistics of the 
incident light is mapped onto the atom-count statistics. Whereas traditional photo-detection theory treats the 
emission centers of photo electrons as distinguishable, here the centers of escaping atoms are condensed and 
thus indistinguishable atoms. From this an enhancement of the photon-number resolution as compared to the 
commonly known counting formula is derived. 



PACS numbers: 42.50.-p, 42.50.Ar, 37.10.Gh 

I. INTRODUCTION 

The quantum theory of photo detection based on the absorp- 
tion of photons and emission of photo electrons represents one 
of the cornerstones of quantum optics. It serves to obtain the 
statistics of emitted photo electrons given the quantum statis- 
tics of the incident optical field. Various approximations, as 
described by this theory, lead eventually to the famous photo- 
counting formula of Mandel (author?) |flll2fl, a quantum ver- 
sion of a previously known semi-classical Poissonian formula 
(author?) Si], 

^„(r,ro) = (:1^^4f^.-''^^^^^^^^ (1) 

Here : : denotes normal operator ordering, tjb is the quantum 
efficiency of the detector, and /(f,fo) is the time-integrated 
light intensity incident on the detector's entrance plane. 

This formula has the well-known limit of a purely Poisso- 
nian photo-electron statistics, if the integration time t = f — fo 
is larger than the coherence time of the incident optical field. 
This integration time represents the response time of the de- 
tector system including the connected electronics to amplify 
the generated photo currents. Thus to observe the statistics of 
the optical field, one must ideally have a fast detector and a 
large coherence time of the optical field under study. 

As already mentioned, to derive the Mandel formula, some 
approximations have to be made. These approximations are 
perfectly justifiable for a solid-state detector device that op- 
erates at not too low temperatures. One crucial assumption 
is the distinguishability of the atoms emitting the observable 
photo electrons. Another approximation is found to consist in 
the perturbative calculus used to obtain joint probabilities of 
photo-electron emissions. Together they lead to the Poisso- 
nian operator form, rather independent of the underlying ab- 
sorption dynamics. 

Consider now a device that operates in a rather different 
regime, that is, it may be cooled down to ultra cold tempera- 
tures in order to behave more quantum than a typical solid- 
state photo detector. For example, let us consider a cloud 
of magnetically trapped Bose-condensed alkaline atoms (au- 
thor?) ijsti floating on the surface of a so-called atom chip 



(author?) i, H. Atoms can now absorb incident photons 
to receive part of the photon momentum, giving them suffi- 
cient kinetic energy to escape from the trap, to subsequently 
be detected, for instance by ionization. 

As such a system is highly degenerate, the emission cen- 
ters of escaping atoms, i.e. the condensed atoms themselves, 
are not distinguishable. Furthermore, a perturbative approach 
to calculate the emission probabilities is hardly appropriate, 
as we may deal with Rabi cycles, where atoms absorb and 
stimulatedly emit photons, thereby returning to the conden- 
sate. Thus, the crucial approximations that led to the Mandel 
formula cannot be applied and thus one may expect a rather 
different counting formula. Such a counting formula connects 
the statistics of escaping atoms to the statistics of the inci- 
dent optical field. For the purpose of unveiling the different 
counting formula we study in the following a model detec- 
tor system using a Bose-condensed gas. Although, it serves 
here merely for demonstrating the differences in the resulting 
counting formula, we suppose that this system also possibly 
may be realizable in current experiments. 

The paper is organized as follows: In Sec. |II]the model of 
the photo detector is introduced. The atom-counting statistics 
is then derived in Sec. |III1 followed by a discussion of its 
features in Sec. II VI Finally, a summary and conclusions are 
given in Sec. |V] 



II. PHOTO DETECTOR MODEL 
A. Mechanism of photo-detection 

Let us assume that a cloud of bosonic atoms is magneti- 
cally trapped in a micro trap implemented on an atom chip 
and being cooled well below the condensation temperature, 
T ^ Tc- We suppose that the trapping potential is highly elon- 
gated into one direction and therefore approximate the system 
as being effectively one dimensional. Furthermore, the atoms 
shall interact nearly resonant with a collimated light field with 
incidence parallel to the elongated trap axis, see Fig. [T] The 
overlap of the transverse mode structure of the light with the 
transverse mode structure of the atomic cloud shall be consid- 
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Figure 1: Outline of a photo detector using trapped Bose-condensed 
atoms on an atom chip. Atoms in the condensate absorb photons 
from the incident light and escape from the trap, being then counted 
in the atom detector. 
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Figure 2: Energy-level scheme of the atoms constituting the photo 
detector. The zero-energy level is taken as the threshold, where the 
continuum of unbound excited states starts. The ground-state level 
at frequency — tOt) is superimposed by the vibrational sidebands of 
frequency V. 



ered as a constant mode-matching parameter, that will deter- 
mine the coupUng strength and thereby the efficiency of the 
detector. 

The resonant electronic transition of the atoms shall be 
formed by two levels, the lower level (ground state) being sub- 
ject to the magnetic trapping, whereas the upper level (excited 
state) being unaffected by the trap. Thus, if all atoms start 
from the ground state, those being excited by absorption of an 
incident photon may leave the trap potential to be detected, 
e.g. by subsequent ionization. Some of the excited atoms, 
however, will be de-excited by stimulated emission and there- 
fore will be subject to the trapping potential again. The elec- 
tronic level scheme including the sidebands generated by the 
trap potential, and the loss of atoms from the trap is depicted 
in Fig. 121 

Given the atoms being Bose-condensed at the initial time 
fo, one may ask for the probability P„(f ,fo) to observe a atoms 
escaping from the trap in the time interval [fo,f], during which 
light is incident on the atomic probe. After this interval the 
detector is reset, i.e. all atoms are cooled again into the con- 
densate mode (i.e. the lowest trap level), to start the counting 
again for an identical time interval T = f — fo. Thus T will be 
the analogue to the usual photo-detector integration time and 



the atom-counting statistics Pa will then be related, in a yet 
unknown way, to the statistics of the incident optical field. 



B. Interaction with the incident optical field 

The Hamilton operator of the complete system including 
the detector and the incident optical field can be decomposed 
into free and interaction part as 



H = Hq + V, 
where the free evolution is governed by 

Hq = Ham + Hat, 



(2) 



(3) 



with Hem being the Hamiltonian of the electromagnetic field 
and //at being the Hamiltonian of the atoms. 

The electric field of the incident optical beam can be written 
as a decomposition of monochromatic modes of wave vector 
k and frequency ck. 



dkEkhe"''' + H.a 



(4) 



where Ej^ denote the rms vacuum fluctuations of the electric- 
field modes, bi^ and are the bosonic photon annihilation 
and creation operators, respectively. As the field polarization 
is selected by the resonant atomic transition, only one type of 
polarization is considered here. Using this expansion the free 
Hamiltonian of the electromagnetic field thus becomes 



H„ 



dkckh'iii^. 



(5) 



The atomic system is described by the bosonic atom-field 
operators 't',(x), where i = g,e denotes the electronic state, 
and that satisfy the commutation relations 



The atomic Hamilton operator reads 



2m 



5,.gU(x) 



(6) 



(7) 



where the trap potential only acts in the electronic ground state 
and reads 



U{x) = 



mv^x^ 

Tl(Oeg, 



(8) 



with V and (Oeg being the trap and electronic-transition fre- 
quency, respectively, and m is the atomic mass. 

For the purpose of diagonalizing the free atomic Hamilto- 
nian, we define the Schrodinger eigen-modes: 



2j2 



2m 



-5,,U{x) 



0,>W =hcoi^„<j)i^„(x). 



(9) 
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Thus (j)f,,„{x) with discrete n = 0,1,.. . and eigen-frequencies 

(Og,„ = nV -0)0, (too = COeg - 7), (10) 

are the harmonic-oscillator eigen-modes corresponding to a 
trapped atom in its electronic ground state. The modes 
(l^e,k{x) °= exp(/A:x) with continuous k G [— cxj^cxj] and eigen- 
frequencies 



Zm 



(11) 



are plane waves corresponding to a free atom in its electronic 
excited state. These modes form two independent orthonor- 
mal sets and obey the standard completeness relations 

^(/.;„(..)f,,„(x') = dix~x'), (12) 

j dk<^l,{x)<^,^k{x) = 5{x~x'). (13) 

Each electronic component of the quantized atomic field 
can now be expanded as 

4>g(x) =^|„f,,„(jc), ^e{x)= dkek^,,k{x), (14) 

n •' 

where the operators g„ and e^., each satisfy again the bosonic 
commutation relations, 



[gn,gl,] = 5,„„, [h^el] = 5{k-k') 



(15) 



Using the expansion (fT4l i the free atomic Hamiltonian d?) re- 
duces to 

^at = 'Y^^O}g^nglgn + / dkUCOe^kelek- (16) 
n ^ 

The interaction between the atoms and the optical field 
reads in dipole approximation 

V = -dK^ j dxY^^{x)E{x)^j{x), (17) 

where d is the transition dipole moment of the atoms and 
is the matching between the transverse modes of electromag- 
netic and atomic field. Using the expansions of electric and 
atomic fields, cf. Eqs (|4]i and ( fT4l i. respectively, this interac- 
tion can be rewritten in optical rotating-wave approximation 
as 



V = / ^1^' / dkhakelignbk^^Jk' -k)+n.2i., (18) 



where the (vacuum) Rabi frequency has been defined as 

dEkKi_ 



(19) 



and the Fourier transforms of the trap modes are defined as 



1 



{k)- ^ 



Axe ''■■^0g,„(jc). 



(20) 



Thus the absorption of a photon of wave vector k transforms 
a ground-state atom in trap level n into an excited-state atom 
with a superposition of wave vectors A:' given by (A:' — k) . 
We may therefore define the annihilation operator of an ex- 
cited wave packet, created from trap level n by absorption of 
a photon of wave vector k: 



ek.n"^ J dk'ei,,^^^^{k' -k). 



(21) 



This relation can be inverted to obtain all operators e<. from 
the set of operators e^.^ „ (n = 0, 1,2, . . .) for a specific wave 
vector 



(22) 



The excited-wavepacket operators satisfy the following com- 
mutation relations 

[4,,n4.„']= j Ax4ln{^)e'^''''^'<^g,n'{x), (23) 

and in particular [e|(.,„,e^,J ~ 1. Using these wave packet op- 
erators, the interaction Hamiltonian can be simplified to 

V = -Y.f dkh^lkeljnbk + H.a. (24) 



C. Single-mode approximation 

If we start from a Bose-condensed gas with all atoms being 
in the lowest trap level n = 0, a cycled electronic transition 
will preferably lead again to the lowest trap level by bosonic 
enhancement. We may thus approximate the ground-state lev- 
els by a single mode, corresponding to the lowest trap level 
and may simplify the interaction Hamiltonian (l24l i to 



V 



dknQ.kelQgQbk + H.ii. 



(25) 



Furthermore, we suppose that the incident optical field is quasi 
monochromatic with wave vector ko, so that only the photon 
operator b/^^has to be kept. Using thus the definitions /j^^ b, 
go g, ^ ^1 and ^ko ^ ^ interaction further sim- 
plifies to 



V = -nQ.e''gb + H.a. 



(26) 



The free atomic Hamiltonian, on the other hand, can be 
written in this single-mode approximation as 



Hat = hCOei^i - TlCOog^g 



(27) 



with the average frequency of the excited wave packet being 
determined by the wave vector of the absorbed photon and the 
momentum spread of the ground-state trap level: 

f y 

We = J dkCOeAtgflik-kot - 4 + (28) 
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We note that the performed single-mode approximation ne- 
glects the dispersion of the excited wave packet, being now 
considered as a propagating plane wave. 

The effective transition frequency between the relevant two 
atomic levels becomes now 



(0', 



2m 



V 

4' 



(29) 



In resonance this transition frequency is compensated for by 
the frequency ck^ of the optical field, from which we obtain 
the resonance condition 



V 
4' 



(30) 



The master equation (|34] | can be written in the form 



(36) 



where according to Eqs (|5]i, ( |26] |, and dZTl l the non-Hermitean 
effective Hamilton operator reads 

i/eff = hckob'^b + hWee''e-hcoog^g (37) 

with a = \a\e'f. 

E. Atomic pseudo spin 



D. Loss mectianism 

The spatio-temporal mode of the excited wave-packet in the 
single-mode approximation is obtained as 



<l>eix,t) 



(31) 



It corresponds to a motion into the direction of the wave vec- 
tor of the previously absorbed photon with the group velocity 
Vg^We/ko given by 



fikp 
2m 



1 + (2t7) 



-2 



(32) 



For the atomic quantum fields we may define the pseudo 
spin operator S with components {S± = 5, ± iSy) 



and 



2 



S- = {s. 



A A 



2 V 2 



1 , 



where the atom-number operator is defined as 



(38) 
(39) 



(40) 



where 77 = ^o^-^^o is the Lamb-Dicke parameter with Sxq = 
^h/{2mv) being the rms position spread of the trap ground 
level. For a typical magnetic trap potential the weak binding 
regime applies, where the Lamb-Dicke parameter is T] ^ 1. 
Thus the conservation of momentum is approximately granted 
and the excited atom compensates for the momentum of the 
absorbed photon. 

If the excited wave packet has moved over a distance ~ Sjcq 
it may no longer be recycled into the electronic ground state 
by stimulated emission of a photon, as the corresponding spa- 
tial overlap will be close to zero. It thus has escaped from 
the trap. The time of flight for this to happen is given by 
T = 5xQ/vg and the corresponding rate for this to happen, 
y— 2Ttx^^, is therefore obtained as 



fikl K . . . 
7=—- 1+ 27] 
m 77 '■ 



nhk^ 
mr\ 



(33) 



The escape of atoms and their subsequent detection may 
be modeled as an incoherent loss of atoms, described by the 
master equation 



(34) 



where the atom loss is modeled by the Lindblad-form part 



A^e^e + eH 



(41) 



The operators S±- satisfy the standard su(2) commutation re- 
lations, = 2Sz and [5;,5±] = ±S±. The total number 
of excitations - excited atoms plus photons - is given by the 
operator 



N = e^e + b'^b = S- + -+ Pb. 

" 2 



(42) 



The atomic system can now be described in the basis of 
Dicke states (author?) H] 



|A,Ae)at = \A-Ae)g<S)\Ae)e, 



(43) 



where A is the total number of atoms, i.e. A|A,Ae)at = 
A|A,Ag)at, and A^ is the number of excited atoms. The corre- 
sponding basis states for the total system can then be written 

as 



\A,N,t^) ^\A,N-n)^t®\n) 



(44) 



where is the total number of excitations, i.e. .;V|A,A^,«) = 
N\A,N,n), and | 

^)em is a photon-number state. 
Using the definitions of the spin operators ( I38H42I I. the ef- 
fective Hamiltonian ( [37] ) can be rewritten as 



dtp[ 



7 epe 



-^{e^e,p} 



(35) 



Heff = i/o - ^'^'Sz - '-^A ~n\£l\ {S+b + S-P) , (45) 
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where for notational simplicity we defined A' = A + iyjl with 
A — cko — (o'gg being the detuning from resonance, and the free 
Hamihonian is identified as 



H'Q = hckoN- 



(46) 



As this free part commutes with the remainder of the effec- 
tive Hamihonian, we may transform into the interaction pic- 
ture with respect to Hq, to obtain the master equation in the 
interaction picture 

d,p = (i/effp - pHeff) + 7epe^ (47) 

where the transformed effective Hamihonian becomes 

Heff = -hA'S,, - ~n\Q.\ {S+b + S-b'') . (48) 

For notational convenience we omitted here any indication of 
being in the interaction picture. 

In the master equation (07]), the last term describes the es- 
cape of an excited-wavepacket atom from the trap. The re- 
sponsible operator e, that annihilates one such excited atom 
from the system, can be written in the basis of the states ( |44] | 



Thus accordingly we may replace the operators in the effective 
Hamihonian ( |48] l to obtain 



hA- 
2 



-ri\^o\ 



H.a. 



where the operator 



J N 1 



(57) 



(58) 



has been introduced. 

For a proper functioning of the detector we assume that the 
number of atoms in the gas is much larger than the maximum 
number of photons of the incident optical field. Thus the oc- 
cupied eigenvalues of the operator M are very large and con- 
sequently we may perform an expansion over a small param- 
eter being proportional to the inverse atom number (author?) 
lUKl^. Thus the interaction part of the effective Hamihonian 
is expanded as 



L^VM 1 



M 



e = VN~^i\A-l,N-l,t7){A,N,n\. (49) 



F. Limit of large number of atoms 



L_ VM 1 



2M 



(59) 



where the expansion parameter (L- — j ) /M is chosen to cancel 
the first-order contribution in dSTT i. obtaining the zero-order 
Hamihonian as 



Let us now consider the action of the operators appearing in 
the effective Hamihonian (l48T l on the basis states (|44] |. Firstly, 
the actions of the operators S+b and S, on a state \A,N,n) are: 



S+b\A,N,n) = ^/{N-n + l){A-N + n)n (50) 
x\A,N,n-l), 



S,\A,N,n) = iN-n--)\A,N,n) 



(51) 



Defining a new spin operator L with its components being 
defined by the actions 



L+\A,N,n) = y^{N-n){n + l)\A,N,n + l), (52) 
L^\A,N,n) = y/{N-n + l)n\A,N,n-l), (53) 



L\A,N,n) - («-^)|A,A^,«), 



(54) 



and obeying the usual angular-momentum commutation rela- 
tions, Eqs ( fSOl l and (fSTT l can be rewritten as 



4ff(f) 



-hA' 



N 



-M- 



(60) 



III. ATOM-COUNTING STATISTICS 

A. Counting statistics 

From the solution of the master equation we need to extract 
the probability for a atoms having escaped in the time inter- 
val [to,t], starting at the initial time fo with a perfect Bose- 
condensed gas with A atoms in the trap ground state. Thus the 
initial state at time to can be written in the form 



P(f0) = |A,0)at(A,0|®Pem(fo), 



(61) 



where j9em(fo) is the density operator of the incident optical 
field at the initial time Jq. The latter may be expanded in 
photon-number states as 



Pera(fo) = I^P„y(fo)|n)em(n'|, 



(62) 



S+b\A,N,n) ^ L^(A-^+L] \A,N,n), (55) 
SM,N,n) = (^l-A-L\\A,N,n). (56) 



so that the complete initial density operator can be written in 
the basis states (l44l i as 



P(fo) = Y,Pn.,n'{tQ)\A,n,n){A,n',n'\. 



(63) 



6 



Given the initial density operator (|63] |. the formal solution 
of the density operator at time t > to can be obtained from the 
master equation in the form 

P(f,fo) = £p„(f,fo), (64) 

a=0 

where Pa{t,to) is the (unnormalized) conditional density oper- 
ator corresponding to the history of the detector system where 
in total a atoms have escaped in the time interval [to,t]. The 
norm of this conditional density operator is the probability for 
this history to occur, which is the desired probability to count 
a atoms escaping from the trap: 

P„(f,fo)=Tr[p„(f,fo)]. (65) 

The conditional density operator P(,(r,fo) is itself a sum of 
all possible histories where a atoms escape in such a way, that 
the jth atom escapes at time tj G [tQ,t], where j = 

Pa{t,to) = [ dt,,... [ t/fip„(f,fo;f„,...,fi). (66) 

JlO Jlo 

I 



It is a non-unitary evolution ^ intermittent by so-called jump 
operators ^ that describe the escape of a single atom from 
the trap. The super operator of the non-unitary evolution is 
defined as 

^(r)p = J7eff(f)pj/Jff(0, (70) 
the effective non-unitary evolution operator being 

J7eff(f)=exp(^-|//eff), (71) 
I 



where the transition amplitudes are 

<;^(r) = ^YiN ~m){A,N,m\Uesit)\A,N,n). (76) 



The norm of the conditional density operator on the rhs is the 
joint probability density for a atoms to escape from the trap at 
times fi, . . . ,r„: 

/7„(f,fo;f„,...,fi) =Tr[p„(f,fo;fa,...,fi)]. (67) 

Thus the required counting statistics is obtained as 

Pait,to)= [ dt„... dtipait,to;ta,---,h)- (68) 



B. Quantum trajectories 

Given the initial density operator ( |63] l, the conditional den- 
sity operator Pa(t_Jo',ta, . . . ,fi) is given by the quantum trajec- 
tory (author?) lillilB 



(69) 



and the escape of an atom is described by the super operator 



/p = Yepel (72) 



Thus the joint probability density ( |67] | can be written using 
Eqs (|69]l-(l721) and the initial state as 



(73) 



(74) 



(75) 

I 

Here we made use of the fact that the effective Hamiltonian 
preserves both the atom number and the total number of ex- 



pa{tMta,---,tl)=J^it-ta)/.yV'{ta-ta-l)/ / [h ~ to) p {to) . 

I 



Pa{t,t0\ta, ■ ■ ■ ,fl) = ^P„y (fo)(«>(f ,fo;f«, • • • ,fl \n')\^{t,to;ta, ■ ■ • ,fl|«)), 

where the (unnormalized) quantum-trajectory states starting with the initial state |A,n,«) are 

|4>(f,fo;f„, . . . ,fi|n)) = /'/^J/eff(f - ta)e...eU^n{tx ~ to)\A,n,n). 
Using the representation of the operator e in the basis states, Eq. ( |49] l, this state vector can be rewritten as 

|4>(f,fo;fa,...,fi|n)) = Uef{{t-ta)\A-a,n-a,na) 

lla "1 
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citations. As the trajectory dTSl l has exactly A — a atoms and 
n — a remaining excitations, in the sum of Eq. ( |73] ) only terms 
with n = n' contribute: 

Pc,{t,to;ta, • • • = ^^'n;'a(r,fo;fa, ■ ■ ■ ,fi |n), (77) 

n 

where f„ = p„,„ (fo) is the initial photon statistics and the prob- 
ability density conditioned on initially « photons is defined as 

Pa{t,to;ta,...A\n) = |l|4>(f,fo;f„,...,fi|«))f . (78) 
Thus the atom-counting statistics can be written as 

Pa{t,to)=T.P«{tMn)Pn, (79) 

u 

where the conditional probability for a atoms to escape in the 
time interval [to,t] given that « photons are present is 

Pa{t,to\n)^ [ dta... dtipa{t,tQ;ta,---,ti\n). (80) 

C. Over damped resonant regime 

The features of the dynamics of the absorption and stimu- 
lated emission of photons and the loss of atoms from the trap, 
given a state with A atoms and excitations, depends on the 
saturation parameter (cf. App. |A]l 



Sa.n — 



A2 + (7/2)2 



(81) 



whereM|A,A?,«) =MA^N\A,N,n) withM^.^v =A-{N-l)/2. 
Given perfect resonance of the incident monochromatic light 
field, A = 0, for Sa.n > 1 Rabi oscillations occur that involve 
cycles of absorption and (stimulated) re-emission of photons 
until an atom is lost from the trap. 

In the opposite over damped case, Sa.n < 1, however, no cy- 
cling transition is observed but photons are absorbed and their 
excitation is removed from the system by an atom leaving the 
trap. In other words, the recoil energy is much larger than the 
effective coupling energy of the atoms with the photon field. 



(fiko 



2m 



>^[2ti\^l\^/M^. 



(82) 



We recall that in our case the Lamb-Dicke parameter Tj > 1. 

Thus the number of absorbed photons equals approximately 
the number of lost atoms, which means that only the transition 
amplitudes 



'i'i'Ut) = Vrin-m){A,n,m\Ue«it)\A,n,n) (83) 

have to be considered, since only they start from initial states 
with all excitation being photonic. From Eq. (iTSl l it becomes 
then clear that among these transition amplitudes we may fur- 
ther consider only those where one photon has been absorbed, 
i.e., 



However, as these are approximations we must re- 
normalize correctly the transition amplitudes to obtain statis- 
tically correct quantum trajectories. Originally the normaliza- 
tion read 



(85) 



meaning that starting from a state \A,N,n) the system even- 
tually will end up with certainty in one of the states \A,N,m) 
with m = 0, . . . ,A^. Taking into account now only the transi- 
tion amplitudes (|84] |. we must use the re-normalized transition 
amplitude 



^'-"i,„(0 - ^Ajt) = >^'^^,„(0/ Ve;,„, (86) 

where the probability for the considered transition is defined 

as 



(87) 



In this way we obtain an atom waiting-time distribution 

WA.nit)=\^'AAt)\\ (88) 

that is properly normalized: 

'dtwA.„it) = l. (89) 



The quantum trajectory in the over damped regime is thus 
simplified to 

|<I>(f,fo;ffl,---,fi|n)> = Ueff{t -~ta)\A-a,n-a,n-a) 

a-l 

X Yl'i'A-k^n-kitk+l-tk), (90) 

k=0 



(91) 



so that the joint probability density ( iTSl l becomes 

Pa{t,to;ta,...,ti\n) = WA-aji-a{t - ta) 



a-l 
k=0 



where 



WA.nit) ^ {A,n,n\U^^ff{t)Ues{t)\A,n,n) 



(92) 



is the probability that no atom leaves the trap within the time 
interval t starting from the state \A,n,n). 

When approximating the transition amplitudes, cf. Eq. 
( l86b . to obtain statistically correct quantum trajectories that 
lead to a normalized atom-count statistics P^, also the prob- 
ability ( |92l ) must be consistently approximated. This can be 
done by using the general relation 



W^,„(r) = l 



dt' WA.nit'), 



(93) 



(84) employing on the rhs the approximation 



8 



Using the above results, the conditional probability for a 
atoms to leave the trap given that n photons are present, Eq. 
dSOl l. becomes now 



Pc,(tMn) = / dta 

J to 

a-l 



' dhWA-a.n-a{t-ta) 



to 



k=Q 



(94) 



This convolution integral is expressed as the inverse Laplace 
transform 



-1 



a-l 



k=0 



(95) 



where t = f — fo. 

In the resonant case (A = 0), low saturation (Sa ^1), and 
assuming numbers of photons much lower than the atom num- 
ber, n ^ A, we obtain the atom waiting-time distribution as 
[cf. Eq. dll)] 

WAAt) - sinh2 (^v/T^) e-Tl"-(«-i)^T^, (96) 

where the saturation parameter ( lAlb has been approximated 
by Sa = 4n^A/(7/2)^. This form shows a behavior quite 
similar to over damped Rabi oscillations of a two-level sys- 
tem interacting with a resonant laser field, with the saturation 
being now dependent on the number of atoms. The individual 
Laplace transforms become then 

m-k.n-k{z) = (97) 
(n — k){n — k + q) {n — k + 2q) 
{n-k+ Toz) {n - k + ToZ + q) {n - k + Tqz + 2q) ' 

where we defined 



(98) 



As we deal with low saturation, S'a <C 1, the introduced pa- 
rameters behave as ^ 1 and w 75a /4. 



Thus the product of transformed waiting-time distributions 
in Eq. (l95T l becomes 

a-l /n\(n+ci\rn+2ci\ 

n \aJ \ a ) \ a ) (QQ\ 
k p ~ ^n+To^~J ^«+<7+To,-:-j ^+2q+T:QZ'^ ' ^ ' 

where the Binomial coefficients are defined by means of the 
Gamma function, 



r(x+i) 



y) T{x^y+\)Y{y+l) 



(100) 



The Laplace transform of the no-escape probability is accord- 
ing to Eq. (|93]l 



lLA—a.n—a\~/ 



^-}VA-a,n-aiz) 



(101) 



which can be written as the sum 

To 



Wa (z) = 



n — a + TqZ 
To(n-a) 



(102) 



(n - fl + Tqz) (« — fl + TqZ + q) 

To (n — fl) (n — fl + q) 
(« - fl + Tqz) {n - a + Toz + q) {n — a + Toz + 2q) 



+ 



Thus the complete conditional probability ( |95] l becomes the 
sum of the three inverse Laplace transforms 



a+i 

n \ fn+q\ fn+lq- 



ora')r.^') 

\u+ll \ u l\ a I , \a+ll\a+l)\ a ) 
I To-' \ (q+T:QZ\ (2q+to:.\ ( Tqz. \ (q+Ta:.\ (2q+ZoZ\ 
\a+l)\a+lJ\ a J Va+1/ V n+1 / V a+l J 



(103) 



Each Laplace transform results as a Meijer G function (au- 
thor?) lul, 



Pa{q,p\n) = {a\)'(" + ^'^]{l-pr 



+ (« + l) 



+ia + iy 



a+l 



"W« + ?^G30,i_^ 

al \ a ■ - ' 



1 ^ + 1 2^+1 

—a q — a+l 2q — a + l 



« + ^30 



j-33 



1-p 



1 q+l 2q+l 
-a q — a 2q — a + I 

1 ^+1 2^+1 
—a q — a 2q — a 



(104) 



where we have introduced the parameter 



IV. DISCUSSION 

(105) xhe effective rate at which an atom leaves the trap is Tq" ' » 
4|i2pA/7 so that the introduced parameter p, cf. Eq. (1105b 
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Figure4: Conditional atom-count statistics P„(^,p|n) for«= lOpho- 
tons, p = 0.9, q = 100 (a) and q = 10 (b). The lower part shows the 
deviation from the corresponding Binomial statistics with quantum 
efficiency r]o = 0.8862 for q = 100 (c) and rjo = 0.7730 for q = 10 
(d). 



Figure 3: Dependence of the scaled quantum efficiency rjo/p on the 
saturation-dependent parameter q for n = I (solid), n = 10 (dashed), 
and n = 20 (dotted) photons. Atom-escape probabilities are p = 0.9 
(a) and p = 0.6 (b). 

can be identified as the probability for an atom to escape from 
the trap by the absorption of a photon. However, there is ad- 
ditionally the large saturation-dependent parameter (7 w 2/5^, 
cf. Eq. ( |98] l. that may change the shape of the conditional 
probability ( 1104b . It is thus not obvious how these two param- 
eters can be merged into a possibly existing single parameter, 
such as the quantum efficiency. 

However, such a quantum efficiency may be defined phe- 
nomenologically, demanding that the average atom count con- 
ditioned on the presence of n photons reads 

fl„ = t7d«. (106) 

Thus, only the fraction r]o of the n photons leads on average to 
a„ escaping atoms. In general, the above relation is not neces- 
sarily linear, as the efficiency may be a function of the photon 
number, revealing a nonlinear relation between incoming pho- 
ton and escaping atom numbers. For our specific case, such 
a phenomenological quantum efficiency additionally depends 
on the atom-escape probability p and the saturation-dependent 
coefficient q, thus we obtain 

riD^riD{q,P,n). (107) 

In Fig. [3] (a) this dependence is shown for an atom-escape 
probability p = 0.9 in dependence of the parameter q for 11= I 
(solid), « — 10 (dashed), and « — 20 (dotted) photons. It can 



be observed that for large values of q all these curves converge 
to the value of the atom-escape probability, i.e. 

lim r]D{q,p,n)= p. (108) 

The same behavior is also observed for a lower escape proba- 
bility p = 0.6, see Fig. [3](b). Thus for sufficiently low satura- 
tion Sa, i.e. for sufficiently large q, a linear regime is attained, 
where the quantum efficiency becomes independent of the in- 
cident photon number « and coincides with the atom-escape 
probability p. 

Does the statistics become then identical to that known 
from the Mandel counting formula? To answer this question, 
we proceed as follows: Given that by use of Eq. ( 1106b we may 
identify a quantum efficiency rjo for each conditional atom- 
count statistics Pa(q,p\n), we may compare it with the corre- 
sponding conditional count-statistics that would correspond to 
the Mandel formula of photo-detection ([T). The latter is given 
by the Binomial statistics, 

Pl.''\m\n) = (^l^r^Ul-riDT-". (109) 

In Fig. |4](a) the atom-count statistics is shown for an incident 
light field with « = 10 photons, with an atom-escape proba- 
bility p — 0.9 and q = 100. From part (c) it can be seen that 
the deviation from the Binomial statistics of corresponding 
quantum efficiency Tjo = 0.8862 is not too large. However, 
this changes if the saturation-dependent coefficient is lowered 
to <7 = 10, see the atom-count statistics shown in Fig. |4](b). 
Now the deviation from a Binomial statistics with correspond- 
ing quantum efficiency rjo ~ 0.7730 becomes rather large, cf. 
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Figure 5: Conditional atom-count statistics Pa{q,p\n) for/; =10 pho- 
tons, p = 0.6, q = 100 (a) and q = 10 (b). The lower part shows the 
deviation from the corresponding Binomial statistics with quantum 
efficiency rjo = 0.5639 for q = 100 (c) and rjo = 0.3968 for q = 10 
(d). 



Fig. |4](d). In fact, the statistics becomes narrower as com- 
pared to the Binomial form. 

Whereas these two cases used a rather large atom-escape 
probability of p = 0.9, for a lower value such as p = 0.6 al- 
ready for ^ = 100 a somewhat larger deviation from the Bi- 
nomial statistics can be observed in Fig. |5](c). However, the 
general trend of larger deviation for smaller values of q is con- 
firmed. 

We may thus conclude that for a given value of p in the 
limit <7 ^ oo we approach the form of a Binomial statistics, 
agreeing with the Mandel formula. This limit corresponds to 
Sa 0, i.e. to extremely low saturation. For such low sat- 
uration the atomic gas absorbs photons one by one, making 
the appearance of bunches of simultaneously escaping atoms 
highly improbable. As in this limit, at a given instant of time, 
no more than one escaped atom may be detected, the indistin- 
guishability of the bosonic atoms can be safely disregarded. 
Thus the Mandel formula is reproduced, as given by its per- 
turbative derivation assuming distinguishable photo-electron 
emission centers. 

For larger saturation, however, it becomes more probable 
that several photons are absorbed simultaneously, generating 
thus bunches of escaping atoms. In this regime the indistin- 
guishabiUty of the bosonic atoms becomes relevant and large 
deviations from the Mandel counting formula are observed. 
They manifest themselves by a reduced fluctuation of atom 
counts as compared to the Binomial form. This noise reduc- 
tion leads to an enhanced photon-number resolution as com- 
pared to a conventional photodetector at equal detection effi- 
ciency, which constitutes an important advantage of a micro- 
trap based detector. 



It should be noted that the detection efficiency rjo deter- 
mines the mapping photon to atom numbers to all orders of 
moments for a conventional photodetector. In our case, how- 
ever, it only connects the mean photon and atom-count num- 
bers. The shape and width of the atom-count statistics addi- 
tionally depends on a second parameter, e.g. the saturation, 
which allows for a difference from the Binomial statistics. 

The present model assumes an ideal Bose gas in the limit of 
zero temperature. However, the inclusion of atomic collisions 
can be consistently made by identifying the lowest trap mode 
as the solution of the Gross-Pitaevskii equation. In Thomas- 
Fermi approximation the only modification will be then the 
increased (decreased) size of the ground-state mode for re- 
pulsive (attractive) scattering, which will modify the time of 
flight of excited atoms to escape from the trap. Thus repulsive 
(attractive) scattering will decrease (increase) the atom-loss 
rate y. 

The treatment of finite temperature (but below the conden- 
sation temperature) requires the departure from the single- 
mode approximation. However, for low saturation as dis- 
cussed here, the inclusion of non-condensed modes will not 
lead to a mixing of condensed and non-condensed atoms by 
photon absorption: Photon absorption will excite both ini- 
tially condensed as well as non-condensed atoms to immedi- 
ately leave the trap. Thus they do not return to the electronic 
ground state, preventing the mixing. The small fraction of 
non-condensed atoms at finite temperature will represent dis- 
tinguishable emission centers that will contribute as a small 
Binomial admixture to the atom-count statistics. With increas- 
ing temperature this admixture will deterioate the enhance- 
ment of photon-number resolution made by the condensed 
atoms, as we expect from the transition beyond the conden- 
sation temperature to a gas of distinguishable atoms. 



V. SUMMARY AND CONCLUSIONS 

In summary we have introduced a model for a photo de- 
tector using a Bose-condensed gas trapped on an atom chip. 
By use of an atomic pseudo-spin approximation based on 
large atom numbers, the photon absorption and atom-escape 
dynamics could be approximated for the over damped case 
where Rabi oscillations are suppressed. A quantum- trajectory 
method then led us to a counting formula involving a sum of 
three Meijer G functions. 

The conditional count statistics has been shown to depend 
on a saturation-dependent coefficient q and an atom-escape 
probability p, to which the quantum efficiency of the detector 
converges in the limit of extremely low saturation, q ^ In 
this limit the differences from a Binomial statistics vanish so 
that our result coincides with the Mandel counting formula. 
However, for realistic values of the saturation substantial de- 
viations can be observed as a reduced atom-count fluctuations. 

An open question is yet the behavior of our model detector 
for saturation 5^ > 1, where full Rabi oscillations start to oc- 
cur. We assume that in this regime we may expect even larger 
and dramatic deviations from the Mandel formula. Then co- 
herent processes must be included, that describe the tempo- 
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rary storage of photon energy in the atomic gas in the absence 
of correspondingly escaping atoms. This is left for a future 
investigation. 
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Using the explicit values for a and j3 we obtain the explicit 
expressions for the above parameters with 5 — i5', 

i\/SA.n sinh(5'/2) 



^1 -5A,„cosh(5'/2) - sinh(572) 

Sa.h — 1 
5A,nSinh'(5'/2) ' 

7T 



5' = Y^l -5a,,„ 



Appendix A: TRANSITION AMPLITUDE 

The transition amplitude can be simplified to 

x(A,«,«-l|e'^l(2|"lv^^-v-.-ric/2)|^^„^„^^ 

where Ma,,! = A — (« — l)/2. The above matrix element can 
be calculated as 

(A,n,«- l|e'(«^=+'^^->)|A,n,«) = V^iz-"'H"'H^ +z)"-\ 
where 

/asin(5/2) 



5cos(5/2)-/j3sin(5/2)' 
52 



z = -- 



a2sin2(5/2)^ 



8 = ^^a^TW 



using a = 2\Q.\Ty^MA „ and )3 = — /7T/2 for the resonant case 
(A = 0). 



where we defined the saturation parameter, 



^A.N 



M^\^MA.^ 

Y\2 



(i) 



(Al) 



In the over damped regime, SA,n < 1, the transition ampli- 
tude becomes thus 



-1 



v/1 -5A,„cosh(572)+sinh(572) 
^ V5sinh(572). 



■i-i 



In the limit Sa.„ ^ 1 this expression reduces in lowest-order 
approximation to 



xsinhi^'^y/ 1 -Sa,„ 



(A2) 
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